Pan‐cancer analyses of bromodomain containing 9 as a novel therapeutic target reveals its diagnostic, prognostic potential and biological mechanism in human tumours

Abstract Background Mutations in one or more genes responsible for encoding subunits within the SWItch/Sucrose Non‐Fermentable (SWI/SNF) chromatin‐remodelling complexes are found in approximately 25% of cancer patients. Bromodomain containing 9 (BRD9) is a more recently identified protein coding gene, which can encode SWI/SNF chromatin‐remodelling complexes subunits. Although initial evaluations of the potential of BRD9‐based targeted therapy have been explored in the clinical application of a small number of cancer types, more detailed study of the diagnostic and prognostic potential, as well as the detailed biological mechanism of BRD9 remains unreported. Methods We used various bioinformatics tools to generate a comprehensive, pan‐cancer analyses of BRD9 expression in multiple disease types described in The Cancer Genome Atlas (TCGA). Experimental validation was conducted in tissue microarrays and cell lines derived from lung and colon cancers. Results Our study revealed that BRD9 exhibited elevated expression in a wide range of tumours. Analysis of survival data and DNA methylation for BRD9 indicated distinct conclusions for multiple tumours. mRNA splicing and molecular binding were involved in the functional mechanism of BRD9. BRD9 may affect cancer progression through different phosphorylation sites or N6‐methyladenosine site modifications. BRD9 could potentially serve as a novel biomarker for diagnosing different cancer types, especially could accurately forecast the prognosis of melanoma patients receiving anti‐programmed cell death 1 immunotherapy. BRD9 has the potential to serve as a therapeutic target, when pairing with etoposide in patients with melanoma. The BRD9/SMARCD1 axis exhibited promising discriminative performance in forecasting the prognosis of patients afflicted with liver hepatocellular carcinoma (LIHC) and mesothelioma. Additionally, this axis appears to potentially influence the immune response in LIHC by regulating the programmed death‐ligand 1 immune checkpoint. For experimental validation, high expression levels of BRD9 were observed in tumour tissue samples from both lung and colon cancer patients. Knocking down BRD9 led to the inhibition of lung and colon cancer development, likely via the Wnt/β‐catenin signalling pathway. Conclusions These pan‐cancer study revealed the diagnostic and prognostic potential, along with the biological mechanism of BRD9 as a novel therapeutic target in human tumours.


INTRODUCTION
SWItch/Sucrose Non-Fermentable (SWI/SNF) constitutes a multi-protein complex tasked with the regulation of gene expression in eukaryotes.It is found to remodel chromosomes mainly by destroying nucleosomes.SWI/SNF is shown to be mutated in a multiplicity of malignant tumours as well as being associated with the survival time of cancer patients. 1SWI/SNF includes three categories: canonical BRG1/BRM-associated factor (cBAF), non-canonical BAF (ncBAF) and polybromo-associated BAF (PBAF). 2 Bromodomain-containing protein 9 (BRD9) is a recently discovered protein coding gene and is a unique member of the ncBAF complex. 3BRD9 includes two domains: a DUF3512 domain and a bromodomain, 4 the latter being the source of the protein's ability to mediate epigenetic modification.Specifically, the bromodomain of BRD9 regulates gene transcription by recruiting the ncBAF complex to the promoter. 5ccording to the results of genomic sequencing, mutations in one or more genes responsible for encoding subunits within the SWI/SNF chromatin-remodelling complexes are found in approximately 25% of cancer patients.In contrast to the extensively studied tumour suppressor genes and oncogenes, which have been the focus of research for decades, the impact of mutations in SWI/SNF genes on cancer remains inadequately understood.In recent years, these genes have gained more rapport, with an increasing number of studies dedicated to exploring the potential therapeutic significance of mutations in genes that encode SWI/SNF subunits.For example, many researchers have focused on their mechanism of promoting cancer and promising targeted therapies. 6There have been some studies showing that BRD9 as a target for anticancer drugs has potential clinical value in some specific cancer types.Specifically, BRD9 bromodomain inhibitors hinder the proliferation of cancer cells while inducing apoptosis in acute myeloid leukaemia, 4,7 rhabdoid tumour 8 and triple-negative breast cancer. 9lthough some existing evidence based on cytological or in vivo experiments supported the association between BRD9 and some specific cancers, there has been no reported comprehensive analysis for BRD9 in pan-cancer.Performing a comprehensive analysis of BRD9 with potential clinical value and evaluation of its association with prognosis and potential molecular mechanisms can offer fresh insights into our understanding of tumourigenesis and the advancement of cancer therapies.In Figure 1, we utilised The Cancer Genome Atlas (TCGA) database to conduct an inaugural pan-cancer analysis of BRD9, encompassing but not limited to gene expression, protein expression, survival status, mutation, protein phosphorylation, DNA methylation, N 6 -methyladenosine (m 6 A) modification, gene set enrichment analysis, immune correlation analysis, etc.We followed up the bioinformatics analysis with experimental validation, employing tissue microarrays from lung and colon cancer patients to detect the expression of BRD9.The biological functions of BRD9 and the signalling pathway involved were also explored in cell lines related to lung cancer and colon cancer.

BRD9 expression profile data analysis
Along with corresponding clinical information, raw RNAseq data counts were acquired from the TCGA portal.The raw data were then processed by converting F I G U R E 1 A workflow of this study.The batch-corrected, normalised expression data and clinical information of The Cancer Genome Atlas (TCGA) pan-cancer and target datasets were retrieved.The comprehensive pan-cancer analysis of BRD9 included gene expression, protein expression, survival status, mutation, protein phosphorylation, DNA methylation, N 6 -methyladenosine (m 6 A) modification, gene set enrichment analysis, receiver operating characteristic (ROC) analysis, immune correlation analysis, a least absolute shrinkage and selection operator (LASSO) regression analysis, etc.For experimental validation, tissue microarrays from lung and colon cancer patients were used to detect the expression of bromodomain containing 9 (BRD9).The biological functions of BRD9 and the signalling pathway involved have also been explored in lung cancer and colorectal cancer cell lines.The experiments contained PCR, cell proliferation, flow cytometry, colony formation, immunofluorescence, co-immunoprecipitation and Western blot.
counts to transcripts per million (TPM) and normalising it using log2(TPM + 1).Tissue samples containing associated clinical information were retained, resulting in a dataset of these samples for subsequent analysis.
BRD9 expression was compared between tumour and their paired adjacent normal tissues across various tumour types present in the TCGA database.The statistical differentiation was assessed by Wilcoxon signed-rank test.The TCGA cancer types lacking corresponding normal sample were excluded from this analysis.
To supplement the data from TCGA database, we also utilised a proteomic database-Clinical Proteomic Tumor Analysis Consortium (CPTAC) 10 to access BRD9's protein expression.The advantage of validation at the protein level is that while most other databases analyse gene expression at the mRNA level, CPTAC can more accurately reflect the physiologic state of the disease.

2.2
The human protein atlas Data on BRD9 expression in cancer cell lines and cancer tissues were obtained from human protein atlas (HPA) (https://www.proteinatlas.org). 13

Survival analysis
Overall survival (OS) and disease-free survival (DFS) of BRD9 across all TCGA cohorts were analysed using 'Survival Analysis' module in GEPIA2, and survival maps and plots were generated. 11Kaplan-Meier curves were generated via the 'survminer' package.For survival map, we examined the impact of BRD9 on survival across various cancer types, employing the Mantel-Cox test for estimation and comparison.In the survival plots, cohorts with high and low BRD9 expression levels were separated based on median values.

Genetic mutation analysis
We entered 'BRD9' in the 'Quick Search' section of the cBioPortal for Cancer Genomics (http://cbioportal.org). 14,15The alteration frequency for all the tumour types represented in the TCGA database was shown using 'Cancer Types Summary' module.Mutation sites within BRD9 were represented in both the schematic diagram of its protein structure and the 3D structure, both of which were accessible within the 'Mutations' module.Kaplan-Meier curves were shown in specific TCGA cancer types, both with and without BRD9 genetic alterations, using the 'Comparison/Survival' module.

m 6 A modification analysis
A heatmap obtained from the 'Exploration-Gene_Corr' panel of TIMER 2.0 was utilised to represent the correlation analysis between BRD9 and nineteen m 6 A regulators in diverse cancer types sourced from the TCGA database.The purity-adjusted Spearman's rank correlation test was performed to generate the partial correlation coefficient (cor) and its corresponding p-value.When m 6 A regulators were mutated, BRD9 expression was shown by the 'Gene_Mutation' panel of TIMER 2.0.The prediction of m 6 A modification sites within BRD9 was carried out using the sequence-based RNA adenosine methylation site predictor (SRAMP) web tool (http:// www.cuilab.cn/sramp/) in full transcript mode with the following settings: analyse RNA secondary structure-no; tissue-generic; and show query sequence as RNA. 16

Experimental validation
Experimental validation was performed in tissue microarrays and cell lines derived from lung and colon cancers, including quantitative polymerase chain reaction (qPCR) analysis, cell proliferation assay, apoptosis analysed by flow cytometry, colony formation assay, immunofluorescence, Western blot, immunochemistry and co-immunoprecipitation (co-IP).The experimental set-up was described in detail in the Supporting Information.

Constructing and validating the prognostic signature in TCGA cohorts
To compare and assess the predictive accuracy of each gene and risk score, timeROC analysis was conducted. 17eature selection was conducted by the least absolute shrinkage and selection operator (LASSO) regression algorithm via 'glmnet' (version 4.1-1) package, with a 10-fold cross-validation approach. 18Univariate Cox proportional hazard regression and log-rank tests were performed.

Spearman's correlation analysis
The infiltration level of immune cells was calculated though the TIMER algorithm using the R package 'immunedeconv'. 19Spearman's correlations between the immune cell infiltration levels and the risk score of the prognostic signature were calculated, and the results were compared using the 'ggstatsplot' package.

The prediction of potential drug compounds
The potential drug compounds were predicted from the Drug Signatures Database (DSigDB) accessed through the Enrichr platform (https://maayanlab.cloud/Enrichr/)according to similar differentially expressed genes. 20
In Figure 2C, compared in different pathological stage, the BRD9 expression showed significant difference (p < .05)only in ovarian serous cystadenocarcinoma (OV) and LIHC.
Immunofluorescence assay identified that BRD9 was highly expressed in most cancer cell lines, such as epidermoid carcinoma cell line A-431, human osteosarcoma cell line U-2 OS, human glioblastoma cell line U-251 MG, human ovarian carcinoma cell line EFO-21 and squamous cell carcinoma cell line SiHa (Figure 2D).Similarly, BRD9 also showed elevated expression in tissue samples from LIHC, LUAD, LUSC, COAD, READ, STAD and renal cancer (Figure 2E).

OS and DFS analysis
We indicated that elevated BRD9 expression was linked to an adverse prognosis, with regards to OS, for adrenocortical carcinoma (ACC), LIHC, mesothelioma (MESO) and sarcoma (SARC) cohorts (Figure 3A).Conversely, low expression of BRD9 was linked to a unfavourable outcome for OS in lower grade glioma (LGG) brain cancers.Figure 3B showed elevated BRD9 expression was associated with an unfavourable prognosis, with regards to DFS, for cancers of ACC and LIHC within the TCGA database.

BRD9 mutation analysis
Figure 4A and Table S2, sourced from cBioPortal, display the genetic alteration status of BRD9 for patient samples.The highest alteration frequency of BRD9 occurs among patients with LUSC, where the primary type of alteration was 'amplification' referring to copy number alteration (CNA).The following six types of cancer with genetic alteration had especially high copy number amplification of BRD9: ACC, lymphoid neoplasm diffuse large B-cell lymphoma (DLBC), KIRC, MESO, pancreatic adenocarcinoma (PAAD) and thyroid carcinoma (THCA).CNA of BRD9 was particularly pronounced in testicular germ cell tumours (TGCT), where copy number deep deletion made up all the genetic alterations to BRD9 in that particular cancer.Figure 4B,C showed the numerous mutation sites of BRD9 among the TCGA cohort, with emphasis on site with the highest alteration frequency (D167G/Y) in bromodomain in the 3D structure of BRD9.Two samples from the TCGA database were found to be mutated at this shared site: one sample from UCEC cases and another sample from LUSC cases.However, the D167G/Y alteration of BRD9 did not show any effect on the clinical survival on the two UCEC and LUSC cases (Figure 4D).

Tumour mutation burden/microsatellite instability analysis associated with BRD9
High tumour mutation burden (TMB) has been shown to be a biomarker for immunotherapeutic response for treatments that inhibit programmed cell death 1/programmed death-ligand 1 (PD-1/PD-L1) ligation in multiple tumours, thereby improving the survival of cancer patients. 21Similarly, microsatellite instability (MSI) has been suggested as a potential marker for immune checkpoint blockade therapy. 22The relation between BRD9 expression and TMB across all tumour types indicated that BRD9 expression was positively related with TMB in ACC, LUAD, MESO, PAAD and skin cutaneous melanoma (SKCM), but negatively related with that in breast invasive carcinoma (BRCA), THCA, thymoma (THYM) and uveal melanoma (UVM) (Figure S1A).BRD9 expression was positively related with MSI in bladder urothelial carcinoma (BLCA), LUAD, LUSC, SARC and STAD, and negatively related with that in DLBC and pheochromocytoma and paraganglioma (PCPG) (Figure S1B).

Phosphorylation site analysis associated with BRD9
Positive phosphorylation sites were illustrated in a schematic diagram of the BRD9 protein, comparing normal tissues to primary tumour tissues (Figure S2A), and variations in BRD9 phosphorylation levels were depicted using box plots for different cancer types (Figure S2B-F).In the LUAD cohort, no significant differences were observed in the phosphorylation levels of BRD9 between normal tissues and primary tumour tissues.
In contrast to normal tissues, the phosphorylation levels of the T103 locus, situated outside the bromodomain domain of BRD9, exhibited an increase in breast cancer and clear cell RCC tissues but a decrease in UCEC and colon cancer tissues.The phosphorylation The relationship between bromodomain containing 9 (BRD9) expression and survival outcome of multiple cancer types in The Cancer Genome Atlas (TCGA) database.According to log10 (hazard ratio [HR]), overall survival (A) and disease-free survival (B) for different patient cohorts were displayed on the survival maps.Kaplan-Meier curves with positive results were shown (p < .05). level of the S568/S588 locus, located outside the DUF3512 domain, also exhibited an increase in breast cancer, colon cancer and ovarian cancer.Additionally, it is noteworthy that in comparison with normal tissues, the S482 locus within the DUF3512 domain of BRD9 displayed a decreased phosphorylation level in breast cancer tissues.

BRD9 DNA methylation analysis
When compared with normal samples, the level of BRD9 methylation differed significantly in 13 cancer types.We observed increased methylation levels of BRD9 in BRCA, KIRC and KIRP (Figure S3A), and conversely reduced methylation levels in BLCA, ESCA, HNSC, LIHC, LUAD, LUSC, PAAD, prostate adenocarcinoma (PRAD), READ and UCEC (Figure S3A).Relative methylation of the eight cancer types with the most significant differences are represented as heatmaps for multiple probes in Figure S3B-I (all p < .0001).

m 6 A modification analysis associated with BRD9
m 6 A, a prevalent RNA modification, plays a role in various biological processes, including tumour progression. 23The expression level of m 6 A regulators are intricately linked to the activity of cancer-related signalling pathways. 24The heatmap in Figure 5A showedthe relationship between BRD9 and 19 m 6 A regulators among various cancer types.Notably, BRD9 expression was positively correlated with 19 m 6 A regulators of BLCA, HNSC, human papillomavirus (HPV)-negative HNSC (HNSC-HPV-), KIRP, LIHC, SKCM, SKCM-metastasis and UCEC.When any of the m 6 A readers (IGF2BP2 and IGF2BP3) or m 6 A writer METTL3 were mutated, BRD9 was significantly highly expressed in cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), BRCA and LUAD (Figure 5B).In contrast, Figure 5C showed that with mutation of any of the m 6 A readers (HNRNPA2B1, YTHDF1, YTHDF2, YTHDF3 and YTHDC2) or m 6 A writers (METTL14 and ZC3H13), expression of BRD9 was F I G U R E 4 Genetic alteration of bromodomain containing 9 (BRD9) in different tumours based on The Cancer Genome Atlas (TCGA) database by using the cBioPortal tool.The alteration frequency with mutation type (A) and mutation site (B) for multiple cancers were shown.The mutation site with the highest alteration frequency (D167G/Y) in the bromodomain of BRD9 (C).The correlation between BRD9 mutation status and overall, disease-specific, disease-free and progression-free survival of uterine corpus endometrial carcinoma (UCEC) and lung squamous cell carcinoma (LUSC) were analysed (D).significantly downregulated in different cancer types.We also computationally identified m 6 A modification sites for the BRD9 gene sequence (Figure 5D).The following four m 6 A sites were identified with very high confidence: 12 936, 20 470, 27 265 and 28 416.With the previous correlations in mind, this would suggest that m 6 A regulators could affect tumour progression in different cancer types by regulating BRD9 expression.

BRD9-related genes set enrichment analysis
BRD9-related gene set enrichment analysis was conducted, and a total of 29 proteins were predicted to interacted with BRD9 (Table S3 and Figure S4A).We then integrated all tumour expression data from TCGA to determine the top 100 genes that exhibited a correlation with BRD9 expression (Table S4).The top four genes were included.Figure S4B indicated that BRD9 expression showed a positive correlation with that of mediator complex subunit 10 (MED10) (R = .57),NOP2/sun RNA methyltransferase 2 (NSUN2) (R = .65),PAP-associated domain containing 7 (PAPD7) (R = .61),and thyroid hormone receptor interactor 13 (TRIP13) (R = .62).In Figure S4C, the heatmap indicated positive correlation of expression between BRD9 and selected targeting genes (MED10, NSUN2, PAPD7 and TRIP13).The intersection of BRD9-binding genes and highly correlated genes showed only one common gene, SMARCD1 (Figure S4D).Subsequently, Kyoto encyclopedia of gene and genomes pathway analysis was conducted on the genes that interacted with BRD9 and exhibited high correlations.The findings showed that 'axon guidance' and 'spliceosome' could potentially be implicated in the influence of BRD9 on tumour progression (Figure S4E).We additionally carried out Gene Ontology (GO) enrichment analysis on the genes that interacted with BRD9 and displayed high correlations.In biological processes analysis, they were found to be very likely be involved in the regulation of ephrin receptor signalling pathway, peptidyl-tyrosine phosphorylation and mRNA splicing via spliceosome (Figure S4F).In cellular components analysis, the genes were found to be mainly located in the nucleus and nucleoplasm.In molecular functions analysis, the genes were shown to have functions such as poly (A) RNA binding, ATP binding and protein binding.

Expression levels of BRD9 in lung cancer and colon cancer from clinical samples
The expression of BRD9 was evaluated by immunohistochemistry. Figure 6A indicated that BRD9 staining in both LUAD and LUSC was notably higher compared to BRD9 staining in normal lung tissues.When comparing the staining of BRD9 in stage IIB lung cancers against normal lung tissues, only staining in LUSC was significantly higher (Figure 6B).Comparison of BRD9 staining between normal lung tissues and stage IIIA lung cancers once again revealed significantly higher staining in both lung cancer subtypes (Figure 6C).
Similar with lung cancer, the staining of BRD9 in colon cancer, including adenocarcinoma, mucinous adenocarcinoma and signet-ring cell carcinoma, were also notably higher than the staining of BRD9 in normal colon tissues (Figure 6D).The patient details of tissue samples and all images taken for immunohistochemistry were shown in Figures S5 and S6.

3.10
The biological functions of BRD9 in lung and colon cancers Comparing BRD9 expression in three different normal lung cell lines (BEAS-2B, MRC-9 and HLF) against each of the four different lung cancer cell lines (A549, ABC-1, LK-2 and EBC-1) revealed significant upregulation of the gene in tumor.(Figure 7A-C).Among the LUAD cell lines, the highest BRD9 expression was observed in ABC-1 cells.Within the LUSC cell lines, LK-2 cells exhibited the highest BRD9 expression.Similarly, comparing the expression of BRD9 between normal colon cells (CRL1459) and five different colon cancer cell lines revealed that BRD9 was upregulated in colon cancer cells (Figure 7D).In COAD cell lines, BRD9 expression was the highest in HT-29 and SW480 cell lines.The four cell lines with the highest expression of BRD9 in each subtype (ABC-1, LK-2, HT-29 and SW480) were then selected for subsequent experiments.
qPCR was performed on the four different cells, which were transfected with siRNAs, resulting in BRD9 knockdown (Figure 7E-H).The cell proliferation was detected, and it found that BRD9 knockdown significantly inhibited the proliferation of ABC-1 and LK-2 cells (Figure 7I,J).BRD9 knockdown also suppressed the proliferation of HT-29 and SW480 cells (Figure 7K,L).According to the results of flow cytometry (Figure 7M), knockdown of BRD9 promoted apoptosis in lung cancer cells (ABC-1 and LK-2) and colon cancer cells (HT-29 and SW480).BRD9 knockdown also resulted in the inhibition of colony formation in lung cancer cells and colon cancer cells (Figure 7N).
To further validate the unique association between BRD9 and SMARCD1, as shown in Figure S4D, we explored the localization of BRD9 and SMARCD1 with immunofluorescence (Figure 7O).We observed a colocalization between BRD9 and SMARCD1, where SMARCD1 was found to be in the nucleus, which was also the predominant location of BRD9.Analysis of the TCGA-LUAD, LUSC and COAD datasets, further support this association BRD9 expression was positively correlated with SMARCD1 expression (Figure 7P).The Western blot results we obtained lend additional credence to this association (Figure 7Q,R), whereby BRD9 knockdown resulted in reduced SMARCD1 expression in four different cell lines.co-IP assay also confirmed that BRD9 can interact with SMARCD1 (Figure 7S,T).Taken together, our experimental validations paint a detailed picture of the interaction between BRD9 and SMARCD1.
Incidentally, it was found that besides affecting SMARCD1 expression, BRD9 knockdown also reduced the expression levels of Met and β-catenin in Wnt/β-catenin signalling pathway (Figure 7U,V).
In summary, BRD9 interacts with SMARCD1.Knockdown of BRD9 can inhibit the progression of lung and colon cancers through Wnt/β-catenin signalling pathway (Figure 7W).

3.11
Potential clinical value of BRD9

BRD9 could potentially function as a novel biomarker for diagnosis across multiple tissue types
The potential diagnostic value of BRD9 in multiple cancers was assessed through receiver operating characteristic (ROC) curve analysis, with the area under the curve (AUC) serving as the evaluation metric (Table 1 and Figure S7).Based on results from Figure 2A, six types of cancer shown significantly differential expression from TCGA pan-cancer analysis including colorectal adenocarcinoma, ESCA, LUAD, KIRC, LIHC and STAD were selected for validation through logistic regression model.Expression profile of BRD9 in each TCGA cohort was treated as training set, and microarray dataset from Gene Expression Omnibus(GEO) corresponding to each cancer type was used as independent validation.It showed the results of ROC analysis across multiple cancer types.Based on the TCGA cohorts, it is noteworthy that the AUC of BRD9 exceed .7 for several cancer types, including COAD/READ (.832), ESCA (.885), LUAD (.785), KIRC (.794), LIHC (.96) and STAD (.935) (Figure S7).Compared with TCGA cohorts, ROC curve analysis based on GEO validation cohorts of corresponding cancer types yielded similar AUC values.Moreover, the data based on confusion matrices under the optimal 'cutoff' from validation set of various cancer types also showed the BRD9-based models can effectively discriminate between adjacent normal samples and cases of malignancy with high accuracy in these six cancer types.Although there were no GEO datasets available for KIRP and LUSC that can be used for valida- tion, higher AUC values of these two cancer types based on the BRD9 expression profile form TCGA cohorts were also identified (Figure S8): KIRP (.769) and LUSC (.839).These results implied that BRD9 may have the potential to be a novel biomarker for diagnosis across multiple tissue types.

3.11.2
BRD9 expression is associated with the abundance of cancer-associated fibroblasts and the expression of immune checkpoint-related genes Immune response is regarded as one of the hallmarks of cancer, 25 and immune infiltration has gradually attracted widespread attention as a prognostic factor. 26Cancerassociated fibroblasts (CAFs), as the most abundant tumour-associated stromal cells in tumour microenvironment (TME), have the capacity to impact cancer progression by regulating the functions of various tumourinfiltrating immune cells. 27e employed EPIC, MCPCOUNTER, XCELL and TIDE algorithms to investigate the potential link between the prevalence of CAFs and BRD9 expression in various cancer types within the TCGA database.We observed a positive correlation between BRD9 expression and the estimated number of CAFs in TCGA tumours of COAD, HNSC and HNSC-HPV−, but conversely, a negative correlation was observed in LGG (Figure S9A).As the important role of CAFs in cancer biology has become clarified, it has gradually received more widespread attention.Targeting CAFs or their secretome can reduce immunosuppression and remodelling of the TME, which provides effective inroads for treatment of cancer.Checkpoint blockade immunotherapy enhances anti-tumour immune response by targeting T cells.Therefore, combined checkpoint blockade immunotherapy and CAF-targeted therapy is a strategy to treat tumours that thrive in fibroblast-rich TMEs. 28SIGLEC15, TIGIT, CD274, HAVCR2, PDCD1, CTLA4, LAG3 and PDCD1LG2 are immune cell-expressed ligands related to immune checkpoint function. 29The correlations between BRD9 expression and the expression of these immune checkpointrelated genes were analysed for 33 cancer types in TCGA database (Figure S9B).We found that the expression of all eight immune checkpoint related genes above was positively correlated with that of BRD9 in PRAD and PCPG.However, a negative correlation was observed between BRD9 expression and the expression of all eight immune checkpoint-related genes in TGCT.

BRD9 has the potential to serve as a promising therapeutic target in patients with melanoma
Preclinical research on BRD9 is constantly advancing.Of all skin-related cancers, melanoma has the highest mortality rate.Unfortunately, immunotherapies or BRAF/MEKbased targeted therapies are only effective for a subset of melanoma patients.
A study has shown that TP-472 as a BRD9/7 inhibitor can inhibit melanoma tumour growth by promoting apoptosis. 30We downloaded GSE179079, the dataset from this study, to conduct in-depth meta-analysis of their results.Compared with control, 342 genes significantly upregulated and 288 genes significantly downregulated F I G U R E 8 Identification of suggested top drug compounds for the similar differentially expressed genes.(A and C) These significant differentially expressed genes were obtained in volcano plot under treatment with TP-472 at either 5 or 10 μM on A375 melanoma cells for 24 h (log 2 |FC| > 1, p < .05).(B and D) The Kyoto encyclopedia of gene and genomes (KEGG) analysis were conducted based on the differentially expressed genes from (A) and (C), respectively.(E) The similar differentially expressed genes was obtained between treatment with TP-472 at 5 and 10 μM on A375 melanoma cells for 24 h.(F) The suggested top drug compounds for the similar differentially expressed genes were predicted.
were found on A375 cells treated with 5 μM TP-472 for 24 h (Figure 8A).In addition to the previously reported findings, BRD9/7 inhibitor is very likely to inhibit the progression of melanoma by regulating p53, PI3K-Akt or adhesionrelated signalling pathways (Figure 8B).Similarly, we identified 756 significantly upregulated genes and 796 significantly downregulated genes on A375 cells treated with 10 μM TP-472 for 24 h (Figure 8C). Figure 8D showed that these significant differentially expressed genes may be implicated in the regulation of p53 or adhesion-related signalling pathways.
As shown in the Venn diagram in Figure 8E, there was an overlap of 412 differentially expressed genes between Figure 8A,C.Enrichr platform was used to identify candidate drugs for these 412 differentially expressed genes.The results from the candidate drugs were generated based on p-value ranking from the DSigDB database.According to the analysis results, etoposide is the prominent drug compound with which a significant number of genes were connected (Figure 8F).In a previous clinical trial, the combination of etoposide and cisplatin lacked sufficient clinical efficacy in the treatment of metastatic melanoma. 31Our findings would indicate that pairing etoposide with the BRD9/7 inhibitor TP-472 would likely improve clinical outcomes in patients with melanoma.

BRD9 could predict the prognosis of melanoma patients undergoing anti-PD-1 immunotherapy
Myeloid-derived suppressor cells (MDSCs) have the capacity to dampen anti-tumour immune responses, facilitate the development of metastases, and led to resistance against immunotherapy. 32PD-1 has been confirmed to be highly expressed in tumour-infiltrating MDSCs. 33A preclinical study has shown that drug molecule targeting MDSCs enhanced the effectiveness of PD-1 blockade in melanoma, thereby augmenting anti-tumour activity. 34IDE algorithm was employed to investigate the potential connection between the predominance of MDSCs and BRD9 expression across different cancer types within the TCGA database.A positive correlation was observed between BRD9 expression and the estimated number of MDSCs for most cancer types including melanoma (Figure 9A).Based on spatial transcriptome data, BRD9 was shown to be highly expressed in tissues of melanoma patients (Figure 9B).
We further explored the potential of combination therapy of anti-BRD9 and anti-PD-1 in melanoma, utilising data from GSE91061.Based on the median value of the prognostic risk grade, the melanoma samples were divided into low-and high-risk groups.The survival of patients in the two distinct risk categories was shown in Figure 9C.The BRD9 relative expression standards were calculated for every patient.The heatmap results revealed heightened BRD9 expression in the high-risk group, whereas it was comparatively lower in the low-risk group.Survival analysis indicated that the low-risk group exhibited a longer OS in comparison to the high-risk group (Figure 9D).The .5-and 1-year ROC curve analyses was conducted to assess the prognostic accuracy of the BRD9 in melanoma patients undergoing anti-PD-1 immunotherapy.The .5-and 1-year AUC values for the risk signature were .742(95% confidence interval [CI]: .658-.826) and .708(95% CI: .621-.795), respectively (Figure 9E).These findings indicated that the risk score, derived from BRD9, has the potential to effectively forecast the prognosis of melanoma patients undergoing anti-PD-1 immunotherapy.

BRD9/SMARCD1 axis may affect the immune response and predict the prognosis of LIHC and MESO patients
As the results shown in Figure S4D, SMARCD1 was the gene most likely to be correlated and interacted with BRD9.Based on the median value of the prognostic risk grade, TCGA-LIHC samples were divided into low-and highrisk groups.Figure 10A illustrated the distribution of risk grades and the survival of the two groups.The relative expression levels of BRD9 and SMARCD1 were calculated for each patient.The heatmap results indicated elevated expression levels of BRD9 and SMARCD1 in the high-risk group, contrasted with lower expression levels observed in the low-risk group.Survival analysis demonstrated that the low-risk group had a longer OS compared to the high-risk group (Figure 10B).To assess the prognostic accuracy of the BRD9/SMARCD1 axis, we conducted ROC curve analyses using the TCGA-LIHC and TCGA-MESO cohorts.For the TCGA-LIHC cohort, the 1-, 3-and 5-year AUC values for the risk signature were .73(95% CI: .672-.789), .65 (95% CI: .584-.716) and .637(95% CI: .557-.717), respectively (Figure 10C).In the TCGA-MESO cohort, the heatmap indicated a positive correlation between the expression levels of BRD9 and SMARCD1 with the risk grades (Figure 10E), and the low-risk group exhibited a longer OS compared to the high-risk group (Figure 10F).For the TCGA-MESO cohort, the 1-, 3-and 5-year AUC values for the BRD9/SMARCD1 axis as prognostic signature were .671(95% CI: .567-.775), .811(95% CI: .733-.889) and .746(95% CI: .654-.838), respectively (Figure 10G).We carried out a comprehensive analysis to examine the relationship between the infiltration level of immune cells and the risk score derived from the prognostic signature of BRD9/SMARCD1 axis, and found a significant positive correlation between the infiltration level of five specific immune cell types (e.g., B cells, CD4+ T cells, neutrophils, macrophage cells and myeloid dendritic cells) and the risk score associated with the prognostic signature of BRD9/SMARCD1 axis in the TCGA-LIHC cohort (Figure 10D).According to the results of the TCGA-MESO cohort (Figure 10H), there was a positive correlation observed between the infiltration level of B cells and the risk score associated with the prognostic signature of BRD9/SMARCD1 axis.However, the infiltration level of neutrophils exhibited a negative correlation with the risk score derived from the prognostic signature of BRD9/SMARCD1 axis.In short, the AUC values showed that the prognostic signature of BRD9/SMARCD1 axis exhibited a favourable discrimination performance in forecasting the prognosis of patients with LIHC and MESO.These findings suggested that the risk score calculated based on BRD9/SMARCD1 axis could accurately forecast the prognosis of LIHC and MESO patients.Our results also suggested that the BRD9/SMARCD1 axis might affect the immune response and could be used to predict the prognosis of LIHC and MESO patients.

DISCUSSION
Recent reports have highlighted the abnormal expression of bromodomain-containing proteins in various cancer types.They have recently become attractive targets for drug discovery.Based on their structural domains, bromodomain-containing proteins are classified into bromodomain and extra-terminal (BET) and non-BET families of bromodomain proteins.The BET family includes bromodomain containing 2, 3, 4 (BRD2, BRD3, BRD4) and bromodomain testis associated (BRDT). 35The results from Oncomine database 36 showed that BRD2, BRD3 and BRD4 are highly expressed in certain cancer types (Figure S10).Although BRD4 has been recognised as a potential target for cancer therapy, 37 none of the BRD4 inhibitors in current clinical trials have received approval from the United States Food and Drug Administration (FDA) for human use due to 'off-target' effects and induced drug resistance. 37These drawbacks highlight a current unmet need to develop novel small molecule inhibitors for translational medical research.
A promising alternative to BRD4 inhibitors is BRD9based targeted therapy, which has made continuous progress since the protein coding gene was recently identified.We summarised the existing small molecule drugs targeting BRD9 in Table 2.For example, Crawford et al. developed GNE-375, a small-molecule inhibitor targeting the BRD9 bromodomain.This inhibitor has demonstrated the ability to suppress the expression of aldehyde dehydrogenase 1 family member A1 (ALDH1A1), which in turn prevents epigenetically defined drug resistance. 38ramer et al. also found that combination therapy of BRD9 inhibitor I-BRD9 and vincristine can be used as an effective therapeutic approach in rhabdoid tumours. 8f particular note among the current repertoire of potential BRD9-based therapies are the two BRD9 degraders: CFT8634 whose investigational new drug (IND) application was submitted 39 ; and the compound FHD-609 whose phase I study is currently recruiting (ClinicalTrials.govIdentifier: NCT04965753). 40lthough initial evaluations of the potential of BRD9based targeted therapy have been explored in their preclinical research for a small number of cancer types (Table S5), more comprehensive study of the diagnostic and prognostic potential, as well as the detailed biological mechanism of BRD9 remain unreported.Our pan-cancer analysis of BRD9 provides a comprehensive framework perspective for further studies and the development of BRD9-based targeted therapies for a variety of different diseases.
Although the non-BET family proteins (e.g., BRD9) lack the extra-terminal domain, they do contain one or more bromodomains.According to data from the HPA database, 47 BRD9 is mainly located in the nucleoplasm (Figure S11).We analysed 19 different isoforms of BRD9 utilising the data sourced from the TCGA database.The expression level (log2(TPM + 1)) of each isoform was shown in violin plots (Figure S12).Isoform usage distribution (0%−100%) was displayed as bar plots (Figure S13).
Figure 2 suggested that the expression levels of BRD9 were generally upregulated in most tumours when compared to their paired normal samples.Recently, some publications have provided experimental evidence that elevated expression of BRD9 has been linked to the progression of tumours.For example, Huang et al. indicated that miR-140-3p directly targeted BRD9 mRNA, leading to the inhibition of its protein translation.This resulted in decreased proliferation in LUSC by downregulating the expression of C-myc. 48Similarly, Dou et al. discovered that BRD9 promoted the progression of LIHC by activating the TUFT1/AKT pathway. 49ima and coworkers analysed BRD9 expression in multiple cancer types 50 but produced some results that were at odds with our own findings utilising the updated TCGA database.In their publication, compared with normal tissues, BRD9 expression was revealed to be significantly upregulated in BLCA, kidney chromophobe (KICH), THCA and UCEC, but there was no significant difference in ESCA.In our latest results, it was observed that BRD9 was higher expressed in ESCA compared to normal tissues.However, when compared to normal tissues, no significant difference was found in BRD9 expression in BLCA, KICH, THCA and UCEC.
The primary biological function of BRD9 is associated with epigenetic modification, which is mediated by its bromodomain.The regulatory mechanism encompasses the recognition of acetylated lysine residues on both histones and non-histone proteins by the bromodomain, and the recruitment of chaperones, contributing to the TA B L E 2 Small molecule drugs targeting bromodomain containing 9 (BRD9).
No. The BRD9 inhibitor was used to identify genes regulated by BRD9 in Kasumi-1 cells involved in oncology and immune response pathways.The BRD9 inhibitor was identified to mimic genetic perturbation of BRD9.LP99 is a potent and selective bromodomain BRD9 and BRD7 inhibitor with greater potency for BRD9 compared with BRD7.40 modulation of gene transcription. 5BRD9 is shown to be recruited to chromatin binding sites.While there is limited information available regarding the precise biological ligands that interact with the BRD9 bromodomain, there have been indications that BRD9 has the capacity to bind to diacetylated H4K5acK8ac and dipropionylated H4K5prK8pr. 51Little research has focused on BRD9's other domain, DUF3512. 5 Currently, our current knowledge is confined to the recognition that the DUF3512 domain is crucial for BRD9's interaction with other subcomplex subunits. 3Besides this gap in knowledge with regards to the structure of BRD9, there is also no evidence to support the biological significance of changes in the phosphorylation level of BRD9 protein in the occurrence and progression of cancer.Our findings in Figure S2 suggested a potential avenue for future research in this area: investigating whether S482 locus within the DUF3512 domain of BRD9 can affect the progression of breast cancer.Following this thread of inquiry could also help reveal the diversity of function of the DUF3512 domain.

Name of drug
Integrating the data of Figures 2 and S3, we found that downregulation of BRD9 DNA methylation levels in tumour samples, resulted in the upregulation of BRD9 expression level in ESCA, HNSC, LIHC, LUAD, LUSC and READ.Conversely, this was found to be positively correlated in KIRC and KIRP.Generally, promoter hypermethylation can inhibit gene transcription. 52In cases of KIRC and KIRP, there was a notable negative correlation observed between BRD9 DNA methylation and gene expression at the probes located within the promoter region.However, very few probes with significant differences (n = 2) were detected in these two cancer types, as shown in Figures S14 and S15.Although DNA methylation generally has the effect of inhibiting gene transcription, some publications in recent years have put forward new viewpoints.For example, research published in Science by Harris et al. showed that DNA methylation can activate gene expression. 53This could be an explanation for the results in KIRC and KIRP, but it would also suggest a more complex mechanism involved in the influence of DNA methylation on BRD9 gene expression, which should be explore further.
Our results in Figure 4D indicated that, in cases of LUSC and UCEC, the D167G/Y alteration of BRD9 did not affect the prognosis of patients.To pursue this line of reasoning further, we explored the landscape of mutation profiles in LUSC and UCEC samples by employing the 'maftools' package.We analysed somatic mutation data obtained from the TCGA database for these two cancer types.Horizontal histogram plots showed known cancer genes to have the highest mutation frequency in patients with LUSC, such as TP53 (78%), TTN (71%), CSMD3 (42%), MUC16 (39%) and RYR2 (37%).The mutation frequency of BRD9, however, is only 1% (Figure S16A).Similarly, the highest mutation frequency in patients with UCEC was PTEN (57%), and the lowest mutation frequency was BRD9 (3%) (Figure S16B).BRD9 has been reported to have the desirable attribute of being stable in blood, 7 lending itself to a role in blood profiling, which is highly desirable due to being noninvasive, easily accessible and cost-effective. 54Figure S7 highlighted the potential of BRD9 as a diagnostic marker across various cancers, opening up the possibility for using blood assays of BRD9 as a novel detection method to improve the accuracy of cancer screening.
At present, there exists a range of approaches for treating cancer in humans, such as chemotherapy, radiotherapy, targeted therapy and surgery.These treatments, whether administered individually or in combination, are effectual but fall just short of the fabled 'cure for cancer'.In recent years, a novel form of anti-cancer therapy has achieved several miracles of 'clinical cure' for advanced cancer cases.This groundbreaking anti-cancer therapy is immunotherapy, with the most remarkable and widely employed form being anti-PD-1/PD-L1 immunotherapy. 55D-1 is primarily found on activated T cells, whereas PD-L1 is reported to be expressed on various types of tumour cells.Blocking the interaction between PD-1 and PD-L1 can enhance T-cell responses and facilitate anti-tumour activity.Agents that block the PD-1/PD-L1 pathway have been documented to exhibit substantial anti-tumour effectiveness in the treatment of cancer patients.56 Some PD-1/PD-L1 immune drugs have even been approved by FDA, such as nivolumab and pembrolizumab (for lung cancer) and Atezolizumab (for breast cancer).57 Despite the promising clinical results, the effectiveness of anti-PD-1/PD-L1 drugs is not consistent.There is a subset of patients who do not experience benefits from anti-PD-1/PD-L1 therapy, which is termed primary resistance.Additionally, some individuals who initially respond to the treatment may later experience relapse, known as acquired resistance.Combination therapy targeting the PD-1/PD-L1 pathway and resistance mechanisms provides the basis for improving the sensitivity of resistant patients.58 Given the restricted therapeutic impact of anti-PD-1/PD-L1 therapy when used in isolation, there is an immediate necessity to investigate effective combination strategies to overcome resistance to anti-PD-1/PD-L1 therapy, which offers valuable insights for clinical implementation.The TCGA has a few limitations.The tumour samples collected for the atlas were from untreated patients.59 In addition, there is no immuno-oncology data.60 We therefore used the GEO dataset to explore the prognostic potential of BRD9 in immuno-oncology.The results from Figure 9 indicated BRD9 could predict the prognosis of melanoma patients undergoing anti-PD-1 immunotherapy, which could lead to improve immunotherapy efficiency in melanoma patients.
Based on the results of Figure 10, we further evaluated whether the prognostic signature of BRD9/SMARCD1 axis can affect the immunotherapy mediated by PD-L1 blockade.Analysis of TCGA-LIHC data revealed that PD-L1 expression was positively related with the expression levels of both BRD9 and SMARCD1 (Figure S17A).In contrast, TCGA-MESO data indicated that PD-L1 expression did not exhibit a significant correlation with the expression levels of either of BRD9 or SMARCD1 (Figure S17B).This result suggested that BRD9/SMARCD1 axis could have a specific involvement in the immune response of LIHC by regulating PD-L1 immune checkpoint.
Patients with liver cancer are frequently diagnosed at an advanced stage, resulting in an unfavourable prognosis.LIHC accounts for over 90% of liver cancer cases with the best options for treatment being chemotherapy and immunotherapy. 61Based on our preclinical research, use of BRD9 as a biomarker could effectively improve the accuracy of LIHC diagnosis.Furthermore, drug development based on BRD9/SMARCD1 axis may improve the effectiveness of immunotherapy mediated by PD-L1 blockade.
Fang et al. found that BRD9 facilitated the progression of hepatocellular carcinoma by activating the Wnt/βcatenin signalling pathway. 62According to our experimental results, BRD9 knockdown could inhibit the progression of lung and colon cancers by the Wnt/β-catenin signalling pathway.As we have shown, knockdown of BRD9 also reduced SMARCD1 expression.It has been reported that SMARCD1 is a subunit of the SWI/SNF complex. 63hese findings have confirmed, quite comprehensively that BRD9 can interact with SMARCD1.
Taken together, our pan-cancer study revealed the diagnostic and prognostic potential, as well as the biological mechanism of BRD9 as a novel therapeutic target in human tumours for the first time.These findings will be valuable for comprehending the biological role of BRD9 in tumourigenesis and the progression of tumours.They can serve as a foundation for future research and the design of pharmacotherapies aimed at targeting BRD9 for therapeutic purposes.

C O N F L I C T O F I N T E R E S T S TAT E M E N T
The authors declare they have no conflicts of interest.

D ATA AVA I L A B I L I T Y S TAT E M E N T
All data generated or analysed during this study are included in this published article and its Supporting Information.

F I G U R E 2
Bromodomain containing 9 (BRD9) expression in multiple tumours and major pathological stages.(A) Compared with paired normal samples, BRD9 expression levels in tumour samples from multiple cancers types were obtained based on The Cancer Genome Atlas (TCGA) data.There were cases where the sample size was less than three or the within-group standard deviation (SD) was 0 (such as sarcoma [SARC], skin cutaneous melanoma [SKCM], thymoma [THYM]), these groups were not included for statistical analysis, but visualisation was performed.No normal samples were available for some cancer types (e.g., adrenocortical carcinoma [ACC], breast invasive carcinoma [BRCA], lymphoid neoplasm diffuse large B-cell lymphoma [DLBC], etc.).(B) At the protein level, BRD9 expression between normal tissue and primary tissue of breast cancer, ovarian cancer, colon cancer, clear cell renal cell carcinoma (RCC), uterine corpus endometrial carcinoma (UCEC) and lung adenocarcinoma (LUAD) were analysed according to data from the Clinical Proteomic Tumour Analysis Consortium (CPTAC) dataset.(C) The BRD9 expression levels were detected according to pathological stage in ovarian serous cystadenocarcinoma (OV) and liver hepatocellular carcinoma (LIHC) (p < .05).Only statistically significant differences were shown.log2(TPM + 1) was calculated for log-scale, where TPM represents transcripts per million.(D) BRD9 expression was detected by immunofluorescence.BRD9 was labelled with green fluorescence.Microtubules were labelled with red fluorescence.(E) Representative images of BRD9 expression in cancer tissue samples were shown.There were two immunohistochemistry images for each cancer types.*p < .05;**p < .01;***p < .001.

F
I G U R E 5 N 6 -methyladenosine (m 6 A) modification analysis associated with bromodomain containing 9 (BRD9).(A) The correlation analysis between BRD9 and 19 m 6 A regulators in different cancer types based on The Cancer Genome Atlas (TCGA) database were shown in the corresponding heatmap.(B) With mutation of any of the 3 m 6 A regulators (IGF2BP2/IGF2BP3/METTL3), BRD9 expression was detected in different cancer types including cervical squamous cell carcinoma and endocervical adenocarcinoma (CESC), breast invasive carcinoma (BRCA) and lung adenocarcinoma (LUAD) (all p < .05).(C) With mutation of any of the 7 m 6 A regulators (METTL14/HNRNPA2B1/YTHDF1/YTHDF2/YTHDF3/YTHDC2/ZC3H13), BRD9 expression was detected in different cancer types (all p < .05).(D) The computational identification of m 6 A modification sites were shown for BRD9 sequence.

F
I G U R E 6 Immunohistochemical expression of bromodomain containing 9 (BRD9) in human lung cancer and colon cancer.(A) Compared with normal lung tissues, the expression of BRD9 were detected by immunohistochemistry (IHC) staining in lung cancer tissue.Especially, the expression of BRD9 were also detected in IIB stage (B) and IIIA stage (C) lung cancer tissues.Scale bar, 200 μm.(D) Compared with normal colon tissues, the expression of BRD9 were detected by IHC staining in colon cancer tissue.Scale bar, 500 μm.LUAD, lung adenocarcinoma; LUSC, lung squamous cell carcinoma; ns, not significance.*p < .05;**p < .01;***p < .001.

F I G U R E 7
The biological functions of bromodomain containing 9 (BRD9) in lung and colon cancers.(A-C) The expression levels of BRD9 were detected by quantitative real-time PCR (qPCR) in four lung cancer cell lines compared with three different normal lung cell lines respectively.(D) The expression levels of BRD9 were detected by qPCR in five colon cancer cell lines compared with the normal colon cell line.In ABC-1 (E), LK-2 (F), HT-29 (G) and SW480 (H) cell lines, the expression level of BRD9 was detected by qPCR after BRD9 knockdown.The cell proliferation assay (cell counting kit-8 [CCK-8]) was measured in ABC-1 (I), LK-2 (J), HT-29 (K) and SW480 (L) cell lines.(M) The apoptosis rate of ABC-1, LK-2, HT-29 and SW480 cell lines were detected by flow cytometry.(N) A colony formation assay was performed in the four cell lines.(O) The co-localization of BRD9 (red) and SMARCD1 (green) was detected by immunofluorescence in the four cell lines.The cells were counterstained with the nuclear probe DAPI (blue).Scale bar, 10 μm.(P) Based on The Cancer Genome Atlas (TCGA)-lung adenocarcinoma (LUAD), lung squamous cell carcinoma (LUSC) and colon adenocarcinoma (COAD) datasets, the correlation between expression levels of BRD9 and SMARCD1 was analyzed.SMARCD1 proteins were detected by Western blot after BRD9 knockdown in lung cancer (Q) and colon cancer (R) cell lines.The interaction between BRD9 and SMARCD1 was detected by co-immunoprecipitation (co-IP) assay in lung cancer (S) and colon cancer (T) cell lines.The expression levels of key proteins in Wnt/β-catenin signalling pathway were detected by Western blot after BRD9 knockdown in lung cancer (U) and colon cancer (V) cell lines.(W) The schematic of signalling pathway was shown.CCK-8, cell counting kit-8.*p < .05,**p < .01,***p < .001compared with the negative control (NC).

TA B L E 1
Logistic regression model for independent validation of bromodomain containing 9 (BRD9) diagnostic ability.

F I G U R E 9
The potential of combination therapy of anti-bromodomain containing 9 (BRD9) and anti-programmed cell death 1 (PD-1) in melanoma.(A) The heatmap showed a correlation analysis between BRD9 expression and the estimated abundance of myeloid-derived suppressor cells based on TIDE algorithm.(B) Based on spatial transcriptome data from SpatialDB database, 12 the expression levels of BRD9 were detected in tissues of melanoma patients.(C) Distribution of risk score, overall survival (OS), OS status and heatmap of the prognostic signatures of BRD9 in the melanoma cohort form GSE91061 dataset.(D) Kaplan-Meier curves of OS for patients with melanoma who received anti-PD-1 immunotherapy based on expression level of BRD9 in GSE91061 dataset.(E) The prognostic signature was shown by the time-dependent receiver operating characteristic (ROC) curve for predicting .5-and1-year survival based on the cohort in GSE91061 dataset.The 'time' in the figure uses 'year' as the unit.p < .05 was considered as statistically significant.CI, confidence interval; HR, hazard ratio.

F I G U R E 1 0
Construction and validation of prognostic signature of bromodomain containing 9 (BRD9)/SMARCD1 axis in The Cancer Genome Atlas (TCGA) cohorts and immune related Spearman's correlations analysis.(A and E) Distribution of risk score, overall survival (OS), and OS status and heatmap of the prognostic signatures of BRD9/SMARCD1 axis in the TCGA-LIHC (A) and TCGA-MESO cohort (E).(B and F) Kaplan-Meier curves of OS for patients with liver hepatocellular carcinoma (LIHC) based on the risk score in the TCGA cohort (A) and with mesothelioma (MESO) based on the risk score in the TCGA cohort (E).(C and G) The prognostic signature was shown by the time-dependent receiver operating characteristic (ROC) curve for predicting 1-, 3-and 5-year survival based on TCGA-LIHC (A) and TCGA-MESO cohort (E).(D and H) Spearman's correlations between the infiltration level of immune cells and the risk score of the prognostic signatures were calculated based on TCGA-LIHC (A) and TCGA-MESO cohort (E).The 'time' in the figure uses 'year' as the unit.p < .05 was considered as statistically significant.CI, confidence interval; HR, hazard ratio.
displayed enhanced selectivity against the BRD7 bromodomain as well as improved pharmacokinetic properties.
BRD9 bromodomain inhibitor, synergises with BET inhibitors to induce apoptosis in triple-negative breast cancer.
of the n-butyl group as an atypical KAc methyl mimetic allowed generation of 31 (GSK6776) as a soluble, permeable and selective BRD7/9 inhibitor from a pyridazinone template.
is a potent and selective bromodomain BRD9 and BRD7 inhibitor with greater potency for BRD9 compared with BRD7.
is a potent and selective BRD9 inhibitor.
drug (IND) application was submitted for the BRD9 degrader.
Yu Chen, Yuanyuan Fu and Youping Deng designed this research.Yu Chen made figures and tables.Yu Chen performed all experiments.Yu Chen wrote and edited this manuscript.Zitong Gao performed external validation of Figures S7 and S8.Yu Chen, Zitong Gao, Isam Mohd-Ibrahim, Hua Yang, Lang Wu, Yuanyuan Fu and Youping Deng reviewed and revised this manuscript.All authors read and approved the final manuscript.A C K N O W L E D G E M E N T S This work was supported by the National Institutes of Health grants (R01CA223490, R01CA230514, -02-PP90) to Youping Deng.This work was also supported by the 2022 American Association for Cancer Research-Merck Cancer Disparities Fellowship (22-40-68-FU) to Yuanyuan Fu.